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A METHOD FOR DETERMINING NODAL ARRIVAL TIMES 


AT THE MOON FROM PRECESSING NEAR-EARTH PARKING ORBITS 

HAVING VARIOUS INCLINATIONS 

By John R. Unangst and Jesse D. Timmons 
Langley Research Center 


SUMMARY 


Some of the effects of precession of near-earth orbits on the timing of 
inplane launchings from such orbits for a nodal encounter with the moon are 
examined. A range of orbit inclinations from l8° to 58° is considered. From 
comparison with the nonprecessing orbit it is found that the orbit precession 
produces more frequent possible launch times and, in most cases, shorter inter- 
vals between possible launch times. 


INTRODUCTION 


If launch of a vehicle from a near-earth parking orbit and rendezvous or 
impact with the moon is desired, it is most efficient to launch in the plane of 
the parking orbit. In order to avoid the penalties associated with making a 
plane change the launch must be timed so that the vehicle will arrive in the 
vicinity of the moon at the same time that the moon arrives at the line of nodes 
of the lunar-orbit plane and the parking-orbit plane. As pointed out in refer- 
ence 1, the precession of near-earth orbits about the earth 1 s spin axis compli- 
cates the timing of such missions and must be accounted for in determining when 
the moon is in the proper position to allow a launch in the plane of the parking 
orbit. 

In reference 1 the effects of parking-orbit precession on the velocity 
required to correct the lunar trajectory for more than anticipated delay times 
in the parking orbit are considered. In the present paper the effects of 
parking-orbit precession on the frequency of arrival of the moon at the line 
of nodes of the lunar -orbit plane and parking-orbit plane are considered. Var- 
ious values of parking-orbit inclination are considered. It is assumed that the 
angular rate of precession of each parking orbit and the angular rate of the 
moon in its orbit are constants. The resulting equations are solved graphically 
and some numerical results are presented. 



SYMBOLS 


In this paper distances are measured in international nautical miles 
(l international nautical mile = I.852OOO kilometers). 

i L inclination of lunar orhit to earth* s equator, deg 

i;LN inclination of lunar orbit at nodal point N, deg 

ig inclination of parking orbit to earth T s equator, deg 

igN inclination of parking orbit at nodal point N, deg 

r mean distance of parking orbit from center of earth, international 

nautical miles 

equatorial radius of earth, 3^^ international nautical miles 
t time, mean solar days 

Xi,Yj,Zi axes forming inertial geocentric equatorial coordinate system 

with Xj directed toward vernal equinox, Zj directed along 
earth* s axis of rotation, and Yj completing usual orthogonal 
triad 




right ascension of ascending node of lunar orbit on earth’s 
equator, deg 


^LP 

a N 

a-S 


right ascension of moon’s position in its orbit, deg 

right ascension of line of nodes formed by lunar-orbit plane 
and parking-orbit plane, deg 

right ascension of ascending node of parking orbit on earth's 
equator, deg 


5 pp declination of moon's position in its orbit, deg 

declination of line of nodes formed by lunar-orbit plane and 
parking- orb it plane, deg 


*^LP angular position of moon in its orbit, measured from ascending 

node of lunar orbit on earth's equator, deg 


P LS 


angle of intersection of parking orbit and lunar orbit at line 
of nodes, deg 


2 



Subscripts : 


o initial conditions 

1,2 ascending and descending node, respectively 

A dot over a quantity indicates the time derivative of that quantity. 


ANALYSIS 


Figure 1 shows a schematic illustration of the geometrical system consid- 
ered. The XpYpZj- axes form a geocentric equatorial inertial coordinate 

system with the Xj,Yj axes lying in the earth's equator. The lunar-orbit 

plane LNB and the parking-orbit plane SNC are shown projected on a segment of 
an earth-centered unit sphere. The line ON is the line of nodes of planes LNB 
and SNC . 


The problem is to determine how frequently the moon arrives at the nodal 
point N. According to the analysis of reference 2 the oblateness of the earth 
causes the plane of a near-earth orbit to precess about the earth's spin axis 
(Zj axis in fig. l) at a rate which depends on the orbital altitude and the 

inclination of the orbit to the earth's equator. In particular, reference 2 
gives the following equation (expressed here in the notation of the present 
paper) for the precession rate: 


“S 


= - 10. 01 — 


3-5 


COS lc 


( 1 ) 


in degrees per day for orbits of small eccentricity (0.05 or less). Thus, N 
will not remain fixed in space but will vary its position, depending on the 
precession of the parking orbit. It should be noted that the moon's orbital 
plane is not fixed in space but is subject to the perturbing effects of the 
earth and the sun. However, these effects have rather long periods and will be 
ignored herein. (For example, the inclination of the lunar orbit to the eclip- 
tic varies periodically between 28 0 35* and l8°19' in 19*6 years.) It should 
also be noted that the assumption of a circular orbit for the moon eliminates 
the effects of the variation of the angular velocity of the moon. This varia- 
tion amounts to about ±2° per day from the mean value per excursion of the moon 
in its orbit. 

It is apparent from figure 1 that, if the times can be found at which the 
right ascension of the moon is equal to the right ascension of N, the problem 
is solved. Figure 2 shows the unit-sphere projection of lunar- and parking- 
orbit planes and the various angles of interest. The approach taken for this 
investigation is to develop equations for the right ascension of N and the moon 
as a function of time, and by simultaneous solution to determine the times at 
which the equations are equal. 
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In figure 3 the spherical triangles LAN and SAN (fig. 3 (a)) and LNP 
(fig. 3 (h)) from figure 2 are shown. By use of the law of sines the following 
relations are obtained: 

From triangle LAN: 


sm il 


cos Ijjx 


sin 6^ sin^ajj - 


( 2 ) 


Similarly, from triangle SAN: 


From triangle LNP: 


sm i c 


cos 


1 SN 


sin 5 W sin^ - Og) 


(5) 


COS It 


COS 1 


LN 


cos 6 


N 




and from triangle SNP: 


cos 1, 


COS 1 


SN 


cos 5 


’N 


(5) 


Substituting equations ( 4 ) and (5) into equations (2) and (3), respectively, and 
reducing yields the following equations: 


tan 6^ - tan i^ sin^a^ - 
tan 6-^ = tan ig sin (a^ - ag^ 




Eliminating tan and solving for tan cljj yields 

tan i^ sin - tan ig sin ag 


tan 




tan i-^ cos cx^ - tan ig cos ag 


If the precession rate of the parking orbit is assumed to be constant, 
equation (7a) can be written 

tan 11 sin <xl - tan ig sin (ag 0 - agt) 
tan «— - — — 


“TT = 


tan i;L cos - tan ig cos(ag^ 0 - agt) 


(6) 


(7a) 


(7b) 
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Equation ( 7 b) gives the right ascension of the nodal point N as a function of 
time. 


It is now necessary to express the right ascension of the moon as a 
function of time. From figure b it can he seen that 


a LP = a L + 


From triangle LMD 


tan - tan T^p cos ip 


( 8 ) 

(9) 


Then 


tan = tan (ap + 

tan aL + tan t]lp cos 
1 - tan aL tan cos ±l 


( 10 ) 


Since the angular velocity of the moon in its orbit is assumed to be constant. 


^LP ^LP, o + ^LP 1 


( 11 ) 


and equation (10) becomes 


tan ccpp 


tan a L + cos i L tan^p^ + ijjpt) 

1 - tan a L cos i L tan^jp^ + T^pt) 


( 12 ) 


Thus equations ( 7 b) and (12) are the desired equations for this analysis. 

In order to determine the times at which equations (7b) and (12) are equal, 
a graphical solution is used. The equations are evaluated individually and 
plotted as a function of time. Points of intersection of the resulting curves 
are the desired solutions. 


As a matter of interest, the angle between the lunar -orbit plane and the 
parking-orbit plane at the nodal point N is also determined. This angle p 

LS 

is a measure of the plane -change maneuver required to place the rendezvous 
vehicle into the moon's orbital plane at N. 

From triangle LSN in figure 3(a) 


cos p =• cos 
LS 


i L cos ig + sin ij, sin is cos (03 - a L ) 


(13) 
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RESULTS AND DISCUSSION 


The results of some numerical calculations using the previously discussed 
equations are shown in figure 5* These results are presented in table I. For 
these calc \ilations, it was assumed that the inclination of the moon's orbit 
remained fixed at 28°. As indicated by the data presented in reference 3> this 
value corresponds to the inclination of the moon's orbit during the calendar 
years 1967 to 1971 to within ±1°. Values of inclination of the parking orbit 
of l8°, 2 6 °, 28°, 30°^ and 38° were chosen. A parking-orbit altitude of 
228 international nautical miles was used (r =■ 3; 672 international nautical 
miles for eq. (ll)), resulting in a range of parking-orbit precession rates of 
6. 3° to 7*6° per day. The angular velocity of the moon is assumed to be 13.19° 
per day. 

Figure 5 presents the variation of the right ascension of the line of 
nodes ajj, the right ascension of the moon app, and the intersection 

angle p for a 60 -day period corresponding to approximately two revolutions 
LS 

of the moon about the earth. It should be noted that two values of Op are 
plotted at any given time. These values of ajj correspond to each end of the 

line of nodes formed by the two orbital planes and are l80° apart. For the cal- 
culations given in table I, it is assumed that, initially, the moon is at its 
ascending node on the earth's equator ^T)pp Q = 0^, that the ascending nodes of 

both the parking orbit and lunar orbit are coincident (o.L,o = ct S,o)^ and that 
these ascending nodes are alined with the vernal equinox = a L,o = 0). 

It should be noted that these initial conditions place the moon in a position 
so that it is alined with the vernal equinox. These initial conditions are 
assumed for illustrative purposes only and, even though this situation may 
never occur, will not affect the generality of the results. 

Examination of figure 5 shows that the nature of the variation of a,p 
with time changes markedly as the inclination of the parking orbit increases 
from values less than ip to values greater than ip. For values of ig < ip 

(figs. 5(a) and (b)), the variation of with time is seen to be cyclic in 

nature and, furthermore, to involve a periodic alteration in the slope of the 
curves. For values of ig > ip, figures, 5(<i) and (e) show that the variation 
of ojj with time always has a negative slope. For the special case of ig =■ ip, 
figure 5(c) shows that Op varies linearly with time. This special case can be 
verified analytically by setting i g = ip in equation 7('b) and differentiating 
with respect to time, remembering that dp = 0. The result of this operation is 

- \ “s C 11 *-) 

It should be noted that since the moon's orbit is considered to be fixed and 
circular, the variation of Opp in figure 5 is the same for all cases. 
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The time increment "between arrivals of the moon at the line of nodes is 
plotted in figure 6 for each arrival and each parking-orbit inclination. Also 
shown is the time increment between arrivals which would occur if the parking 
orbit were not precessing - 0^ . Figure 6 shows that the precession of the 

parking orbit results, for a majority of the calculations, in reduced time inter- 
vals between arrivals of the moon at the line of nodes when compared with the 
nonprecessing case. For the cases considered herein, a majority show a time 
increment between arrivals of 10 to 11 days. This increment compares with 
I 5.9 days for the nonprecessing case. The longest and the shortest time incre- 
ments are seen to be 21.5 days and 0.9 day, respectively. It will also be noted 
that the number of arrivals during the 60 -day period considered increases from 
four for the nonprecessing case to a maximum of seven for ig - 30°. Thus, for 
the particular examples considered and the initial conditions assumed, the pre- 
cession of the parking orbit results in more frequent possible launch times and, 
in most cases, shorter time intervals between the possible launch times. Fig- 
ure 6 also indicates that parking-orbit inclinations greater than the lunar -orbit 
inclination are more favorable from the standpoint of launch frequency. In par- 
ticular, for a parking-orbit inclination of 30 °* seven possible launch times 
occur during the 60-day period considered. For parking-orbit inclinations 
of 28 ° (i g - i L ) and 38 °, six possible launch times occur. This compares with 

five possible launch times for parking-orbit inclinations of 18 ° and 26 °. 


CONCLUDING REMARKS 


In this investigation the effects of the precession of near-earth orbits 
due to the oblateness of the earth on the timing of launchings from such an 
orbit for a nodal arrival at the moon have been examined. The results of some 
representative calculations performed with the equations derived indicate that 
the precession of the parking orbit results in generally shorter time intervals 
between possible launchings and more frequent possible launch times during the 
60-day period considered. The results also indicated that parking orbits with 
inclinations greater than the lunar-orbit inclination are more desirable since 
these orbits result in more frequent possible launch times and, in most cases, 
shorter time intervals between possible launchings than do parking orbits with 
inclinations less than that of the lunar orbit. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 28, 196 k. 
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Figure Continued. 











Time, t, days 
(e) i s = J8°- 


Figure 5* ~ Concluded. 
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